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HYMO  1987  0005  REPORT 


1.  Introduction 

This  study  relates  to  the  further  development  of  an  operational  model  for 
ungauged  catchment  forecasting. 

The  work  was  begun  under  DAJA  45-83-C-0029  and  continues  under  the  present 
contract.  Figure  1.1  shows  the  overall  research  design  as  viewed  by  the 
authors  In  January  1985.  This  research  design  has  been  modified  to 
that  shown  In  figure  1.2. 

Prior  to  this  report,  validation  and  model  configuration  of  Model  11  (figure 
1.2)  had  been  confined  to  the  following  conditions: 

(I)  single  sub-basin  areas 

(II)  absence  of  channel  routing  therefore 

(ill)  single  time  base  for  precipitation  events  (l.e.  no  moving  storm 
conditions  had  been  evaluated) . 

This  report  details  investigations  that  have  been  undertaken  to  extend  the 
validation  and,  therefore,  application  of  Model  11  into  these  areas. 
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Figure  1.2  :  Overall  Research  Design 
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4. 


Figure  1.3  :  Research  Design  Detail 
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Within  the  context  of  the  background  described  in  Section  1  and  Figures  1.1 
and  1.2,  the  following  objectives  were  specified  for  work  reported  here: 

(1)  An  evaluation  of  the  current  potential  for  the  inclusion  of  land-use 
into  Model  11. 

(ii)  An  evaluation  of  moving  storm  effects  for  subsequent  identification 
of  rainfall/ radar  input  needs  for  given  circumstances. 

(ill)  An  evaluation  of  the  method  Model  11  currently  adopts  for  estimation 
of  the  in  and  out  of  bank  roughness 

(lv)  The  specification  of  the  broad  research  programme  for  the  next  six 
months  in  relation  to  the  conceptual  and  operational  needs  shown  in 
Figure  1.2 

The  basic  model  structure  of  Model  11  was  established  on  the  basis  that  the 
replacement  of  the  Curve  Number  routine  by  the  finite  difference 
infiltration  scheme  offered  higher  resolution  of  model  output  with  little 
change  in  data  input  needs.  This  has  been  broadly  vindicated  by  the 
validation  tests  undertaken  by  Anderson  and  Howes  (DAJA45-83-C-0029,  1986)  at 
the  aub-catchment  level. 

The  development  of  Model  111  (Figure  1.2)  now  demands  that  other  components  of 
the  model  are  examined  and,  if  necessary,  replaced  and  upgraded  in  a  similar 
manner,  consistent  with  operational  demands  and  constraints.  Figure  2.1 
illustrates,  in  general  terms,  a  selection  of  alternative  approaches  within 
each  of  the  model  components.  It  is  clear  from  that  figure  that  we  are 
currently  seeking  to  upgrade  other  modules,  where  appropriate,  in  the  current 
contract. 

Consistency  of  module  compatibility  in  terms  of  parameterisation,  computer 
time  and  output  resolution  is  thus  of  major  concern  in  the  development  of 


Figure  2.1  Current  capability  In  each  model  component 
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Model  111.  (Given  Che  basis  of  the  development  of  Model  11  this  was  much  less 
of  a  factor  of  course) . 


In  particular,  we  see  the  need  to  examine  the  treatment  of  channel  roughness 
currently  employed  In  the  model.  Significant  Improvements  we  believe  are 
possible  In  this  module,  and  this  report  details  Initial  approaches  we  are 
making. 

A  German  watershed  has  been  established  as  a  validation  base  for  much  of  the 
work  currently  being  undertaken,  and  collaboration  with  USDA,  Beltsvllle, 
Washington,  has  resulted  In  a  potential  data  base  for  Oklahoma  being 
Identified  for  detailed  moving  storm/ channel  routing  validation. 
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Evaluation  of  current  potential  for  inclusion  of  land  use 


The  current  capability  of  HYMO  2  is  limited  with  respect  to  vegetation/ land 
use  effects  as  follows: 

(i)  Limited  to  non-forested  areas 

(ii)  Temporally  fixed  hydraulic  conductivity  conditions 

(e.g.  no  soil  crusting  development) 

Given  the  present  input  parameters  (Table  3.1)  it  is  possible  to  envisage 
tackling  this  problem.  To  retain  this  limited  data  input,  and  simultaneously 
incorporate  land  use  specification,  three  options  are  available: 

(i)  To  modify  the  hydraulic  conductivity  for  given  land  use 
conditions 

(ii)  To  modify  the  detention  capacity  similarly,  or 

(iii)  To  modify  the  suction-moisture  curves  in  accordance  with  land  use. 

In  the  context  of  (Hi),  this  facility  already  exists  through  the  inclusion  of 
organic  matter  in  the  Brakenseik  and  Rawls  scheme  we  have  adapted  for  suction 
moisture  curve  generation  (see  report  DAJA-37-81-C-0221 ) . 

No  real  equilibrium  for  (i)  and  (ii)  is  available  at  this  stage,  and  so  we 
have,  in  this  report,  explored  the  sensitivity  of  the  output  hydrograph  to 
change  in  organic  matter. 

2 

A  hypothetical  catchment,  area  24  km  ,  stream  length  11.8  km,  with  no 
subcatchment  subdivision,  was  established  to  test  for  sensitivity  to  organic 
content  for  each  of  the  three  soil  classifications,  clay,  silt  and  sand.  The 
catchment  was  assuned  to  be  spatially  homogeneous,  with  organic  content 
varying  from  0  to  10Z  for  each  class,  in  five  increments,  giving  a  total  of  18 
hypothetical  situations. 


The  other  variables  required  by  the  infiltration  algorithm  were  kept  constant 
wherever  possible  but  some,  like  the  saturated  moisture  content  and  saturated 
conductivity,  vary  with  soil  type.  The  initial  soil  moisture  values, 
restricted  by  the  algorithm  to  be  less  than  the  saturated  value,  were  set  at 
0.04  less  than  the  smallest  saturated  value  in  the  profile  for  ten  cells  in 
all  three  layers.  Hie  last  value  in  the  soil  moisture  characteristic  curves 
is  the  saturated  moisture  value,  and  this  was  matched  with  a  corresponding 
suction  value  of  -0.4  m,  for  all  layers,  and  soil  types.  This  value  was  set 
arbitrarily  but  is  within  the  range  of  bubbling  pressures  (or  soil  air  entry 
points),  determined  by  Rawls  et  al.  (1982).  Detention  capacity  and 
evaporation  were  assumed  to  be  zero  for  all  soil  types  and  storms.  Three 
storms  were  used  with  emulative  rainfall  totals  ranging  from  43  mm  to  205  mm; 
all  were  given  an  identical  duration  of  9  hours. 

Variability  in  hydraulic  conductivity,  initial  moisture  content,  moisture 
content  on  the  suction  moisture  curve  and  detention  capacity  were  included. 

Thus,  for  each  organic  content  value,  within  each  soil  type  and  for  each  of 
the  three  storms,  the  model  was  run  twenty  times;  some  360  times  in  total. 
From  each  run,  the  runoff  volute,  peak  discharge  and  time  to  peak  f~om 
commencement  of  storm,  were  taken  as  characteristics  of  the  hydrograph.  The 
twenty  results  were  used  to  calculate  the  coefficient  of  variation  for  each  of 
the  soil/hydrograph  characteristics/ organic  matter  percentages  shown  in  Table 
3.2.  From  this  table,  the  relative  insensitivity  of  the  hydrograph 
characteristics  to  organic  matter  change  can  be  observed. 

Thus,  land  use  inclusion  by  this  variable  has  a  much  lesser  effect  than,  say, 
the  hydrograph  generation  routine  on  the  output  hydrograph.  To  facilitate  a 
full  land-use  incorporation  into  the  current  approach  would  thus  necessitate 
the  total  restructuring  of  the  current  infiltration  algorithm  currently 
adopted  to  Include  a  physically  based  component  modelling  vegetation  and  soil 
water  interaction  explicitly.  In  the  context  of  mode  component  comparability 
(Figure  2.1),  we  feel  that  this  cannot  be  justified.  An  additional  element 
here  is  the  substantially  Increased  data  that  would  be  needed  to  drive  such  a 
scheme. 
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Table  3 . 1  Data  requirements  of  the  soil  water  model 


Soil  Profile  Hydrologic  Characteristics 
% 

For  each  layer: 

soil  water  content  at  saturation 
saturated  hydraulic  conductivity 
-  suction  moisture  curve  (a  maximum  of  20  observations) 


For  each  cell: 

initial  soil  water  content 


Soil  Profile  Dimensions 


total  number  of  cells  in  column 
number  of  cells  in  layer  1 
number  of  cells  in  layer  2 
thickness  of  each  cell 


Surface  Conditions 


detention  capacity 

maximum  evaporation  during  the  day 


Precipitation 

-  rainfall  data  time  increment 

rainfall  data  for  each  time  increment 
rainfall  start  time 
rainfall  stop  time 


Program  Controls 

-  iteration  time  for  simulation 
simulation  start  time 

-  simulation  stop  time 

number  of  profiles  for  the  catchment  area 


Note:  no  historical  flow  data  is  required. 


Table  3.2  Coefficients  of  variation  based  on  20  runs  for  each  organic  matter/ storm  combination  shown 
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4.  Experimental  catchments  and  related  data  sources 


For  the  purposes  of  model  development  and  for  this  report  the  identification 
of  moving  storm  effects,  test  catchments  in  one  area  have  been  established, 
and  a  further  set  of  catchments  are  to  be  incorporated  into  the  study  within 
the  next  six  months. 

4.1  West  Germany  :  River  Haune,  Fulda 

The  initial  travelling  storm  effects  reported  later  were  developed  on  the 
Haune  catchment  in  West  Germany  (grid  ref.  350  561).  This  is  a  catchment  of 
148  km.  Figures  4.1  -  4.4,  for  which  soil  classification  data,  precipitation 
information,  outflow  hydrographs  and  rating  curves  are  available. 

2 

Radar  precipitation  input,  with  a  resolution  of  5  km  ,  has  been  anticipated. 
Figure  4.2,  as  a  future  development  to  replace  raingauge  data  in  suitable 
areas.  The  effects  of  radar  uncertainty  will  be  Included  stochastically  at  a 
later  date. 

Subdivision  of  the  Haune  catchment,  carried  out  manually,  is  restricted  to 
seven  catchments.  This  limitation  is  Imposed  by  the  need  to  maintain  a 
consistent  level  of  resolution  throughout  the  model  (Figure  2.1),  and  by 
central  processor,  computing  restrictions.  Increasing  the  resolution  of 
subcatchments  in  the  Haune  catchment,  would  extend  the  number  of  catchments 
from  7  to  approximately  25.  Sensitivity  to  subcatchment  resolution  will  be 
examined  in  future  work. 

4.2  West  Germany  r  River  Fulda 

It  is  hoped  to  expand  the  study  site  in  West  Germany  to  include  the  greater 
catchment  of  the  Fulda  river,  of  which  the  Haune  is  a  tributary.  The 
approximate  area  of  this  catchment  is  2500  sq.km.,  about  an  order  of  magnitude 
greater  than  that  of  the  Haune' s  catchment. 


Location  and  selected  routing  reaches 
of  the  River  Haune  used  in  the  moving 
stuim  study  -  section  5. 


4.2  :  Subcatchment  boundaries  < 
radar  grid  (5  km)  used  fi 
inputs  in  the  moving  sto: 


15. 


2 


Kilom#tr«s 


SOILS  DESCRIPTION 


16 


4.3  OSA  :  Oklahoma 


A  further  study  site  is  being  explored  In  Oklahoma,  USA.  This  site  Includes 
some  175  ralngauges  In  2500  sq.km.,  with  which  It  Is  hoped  to  Investigate  the 
modelling  of  severe  convective  storms,  traversing  the  catchment. 

Collaboration  with  U.S.D.A.  at  Beltsvllle,  Washington,  Is  continuing  with 
respect  to  data  acquisition. 

4.4  Hydraulic  research  ;  Hydraulics  Research  Ltd.,  England 

Data  from  hardware  models  and  study  sites  throughout  the  U.K.  is  at  present 
being  Investigated,  In  order  to  Include  out  of  bank  roughness  effects  In  the 
model.  This  Is  discussed  further  In  section  6. 
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The  River  Haune  catchment  was  selected  for  running  HYMO  2  with  precipitation 
input  assigned  to  correspond  to  5  km  grid  squares  (see  Figure  4.2)  to 
replicate  moving  storms,  at  variable  speeds  from  four  orthogonal  directions. 
The  basic  purpose  of  this  design  is  to  be  able  to  evaluate  the  space/ time 
precipitation  input  data  needs  for  HYMO  2  in  the  context  of  hydrograph 
predictions  for  particular  applications. 


The  exact  experimental  design  used  was  as  follows: 

3  storms  -  V 2”  in  24  hours 

1"  in  6  hours 

2"  in  30  minutes 


2  storm  speeds  -  5  km  hr 

10  km  hr”1 


4  directions  -  N,  S,  E,  W. 

Hence,  24  experimental  frames  were  undertaken  in  all.  The  operational  rule  of 
particular  note  in  the  design  adopted  relates  to  the  fact  that  as  soon  as 
precipitation  occurred  in  any  part  of  a  subcatchment,  defined  by  the 
precipitation  grid  square,  then  the  appropriate  precipitation  was  apportioned 
to  the  whole  subcatchment. 

The  results  presented  in  this  report  relate  principally  to  the  1/2"  in  24  hour 
and  1"  in  6  hour  storms.  Table  5.1  illustrates  all  eight  experimental  frames 
for  the  1/2"  storm,  showing  the  time  to  peak  and  peak  discharge  values.  Table 
5.2  details  the  same  parameter  for  the  1"  storm.  The  principal  points  to  note 
from  the  results  are: 

(i)  For  the  1/2"  in  24  hour  storm  (being  typical  of  a  European  frontal 
precipitation  event)  there  is  essentially  no  difference  in  either 
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Hydrograph  results  from  running  storm 
(1"  In  6  hours) 
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time  to  peak  or  peak  discharge  for  any  of  the  conditions  prescribed 
In  terms  of  both  speed  and  storm  direction  (see  Figure  5.1). 

2 

We  may  reasonably  conclude  therefore  that  a  catchment  area  of  145  km 
may  be  regarded  as  the  minimum  for  which  high  resolution  (5  km  grid 
square)  precipitation  data  Is  required  for  low  Intensity  storms. 

For  the  1"  in  6  hour  event,  by  contrast,  certain  differences  are 

apparent  in  the  context  of  both  time  to  peak  and  peak  discharge. 

However,  these  differences  are  not  large  -  see  Figures  5.2  and  5.3. 

Amongst  the  largest  difference  Is  to  be  found  in  the  10  km  hr-*, 

northerly  moving  storm  compared  with  the  5  km  hr  *  southerly  moving 

3  -1 

storm.  In  this  case,  the  peak  discharges  are  respectively  123  m  S  , 
3  —1 

and  114.19  m  S  .  The  respective  time  to  peak  differs  by  some  30 
minutes  -  see  Figure  5.2. 

2 

For  higher  intense  storms  at  145  km  area  then  we  may  expect  maximum 
differences  of  perhaps  no  greater  than  5 Z  in  outflow  hydrograph 
parameters  attributable  to  moving  storm  conditions. 
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Figure  5.1  ;  Hydrographs  for  the  in  24  hour  storm 
(5  km  hr-1)  generated  by  storm  movements 
in  four  orthogonal  directions,  Note 
similarity  of  hydrograph  responses 
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Figure  5.3 

:  Hydrographs  for  the  1"  in 
shovring  for  the  5  km  hr-* 
i)  northerly  storm  track 
ii)  southerly  storm  track 
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6.  Revision  of  out-of-bank  roughness  and  stage-discharge  calculations 

HYMO  2  accepts  that  friction  varies  with  stage  and  Includes  a  function  to 
reduce  Manning's  n  with  rising  stage. 

n'  -  n  -  0.0025  Area  6.1 

Wetted  perimeter 

There  seems  no  physical  basis  for  this  equation  and  with  a  large  hydraulic 
radius,  n'  may  assume  a  negative  value,  resulting  in  a  negative  discharge. 

In  the  flow-rating  curves  generated  by  Manning's  equation  (in  cases  where  they 
are  unknown) ,  there  is  no  provision  for  any  interaction  between  computations 
undertaken  in  the  specified  channel  sections;  l.e.  the  channel  cross 
sectional  sub-areas  are  assumed  independent  in  terms  of  the  hydraulic 
calculations. 

There  are  thus  two  areas  of  potential  research  here: 

(1)  the  assignment  of  roughness  to  channel  sub-areas,  and 

(11)  the  revision  of  the  Manning  equation  for  discharge  calculations  in 

compound  (out-of-bank)  channel  conditions 

Evidence  is  available  to  show  that  significant  errors  can  occur  in  discharge 
estimation  in  compound  channels  using  the  Manning  equation. —  As  we  have  noted 
above,  this  is  due  to  the  in-channel  and  out-of-bank  interaction.  Certain 
workers  (e.g.  Chow)  have  suggested  that  a  redefinition  of  the  area  and  wetted 
perimeter  terms  may  be  undertaken  to  provide  an  improved  discharge  estimate 
based  on  the  Manning  formula.  Figure  6.1  illustrates  the  condition.  Figure 
6.1a  shows  the  full  interaction  between  ln-channel  and  out-of-channel  flow. 
Figure  6.1b  defines  the  channel  morphometry  terms  that  are  then  used  to  define 
the  area  (A)  and  wetted  perimeter  (P)  for  both  the  flood  plain  (Fp)  and  main 
channel  (Me). 
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From  Manning's  equation: 


Qt  -  (At  Sf  )  .  (Pi  n^-1  6.2 

where  is  the  flow  through  each  sub  area  A^.  Table  6.1  provides  four 
different  methods  used  to  determine  A^  and  ,  based  upon  momentum  transfer 
between  sub  areas  (Figure  6.1a).  Knight  (1984)  has  shown  that  whilst  certain 
of  these  methods  provide  differences  in  discharge  estimates  by  up  to  as  much 
as  +25%  of  measured  flow,  method  3  was  significantly  more  accurate  in  the 
fluae  trials  undertaken. 

We  will  be  evaluating  the  effect  of  including  these  alternative  procedures  in 
the  estimation  of  the  rating  curves  for  HYMO  2  in  the  next  six  month  period. 

Evidence  from  flume  studies  conducted  at  Hydraulics  Research  shows  that  a 
pronounced  hysteresis  can  occur  in  the  stage-discharge  relationship  for 
compound  channels.  Potentially,  this  can  be  thought  to  have  a  major  effect 
upon  flood  Inundation  levels,  and  is  currently  not  Included  in  the  HYMO  2 
computation  procedure.  When  the  full  effects  of  equation  6.2  and  Table  6.1 
have  been  evaluated,  this  aspect  will  be  examined. 
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7.  Discussion 

The  following  progress  has  been  made  In  this  reporting  period: 

‘\ 

(1)  A  West  German  catchment , (the  River  Haune,  Figure- 4.1)  has  been 
established  as  a  validation  base  for  HYMO  2. 

(ii)  On  this  catchment,  an  initial  study  has  been  completed  of  the 
relative  effects  of  many  storms  on  resultant  catchment  outflow 
hydrographs  '{section  5).  This  study  is  an  essential  element  of  the 
model/ data  development  programme.  We  aim  to  specify  the 
precipitation  input  needs  (in  time  and  space)  for  selected 
applications  of  HYMO  2  within  the  next  twelve  months. 

(ill)  An  initial  study  of  an  approach  to  land  use  incorporation  into  the 
infiltration  algorithm  using  organic  matter  change  has  revealed  a 
relative  insensitivity  to  this  parameter  (Table  3.1).  To  apply  HYMO 
2  to  forested  conditions  may  necessitate  further  exploration  -  most 
probably  in  the  modification  of  the  hydraulic  conductivity  under  such 
conditions.  It  is  to  be  noted  that  validation  of  our  scheme  so  far 
has  been  confined  to  non-forested  areas  (see  report 
DAJA-45-83-C-0029) .  We  plan  to  rectify  this  within  the  next  twelve 
months . 

(iv)  The  current  treatment  in  HYMO  2  of  rating-curve  generation  using  the 
standard  Manning  equation  has  been  examined.  It  is  considered  that 
the  current  approach  (equation  6.1)  may  lead  to  errors  in  both  peak 
flow  and  flood  inundation  estimates.  An  examination  of  alternative 
methods  of  estimating  area  and  wetted  perimeters  in  compound  channels 
(s«e  Table  6.1)  will  be  examined  in  the  context  of  the  Fulda 

catchment  in  the  next  reporting  period.^ 

\ 
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